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We study the spectral properties of a class of random matrices where the matrix elements depend 
exponentially on the distance between uniformly and randomly distributed points. This model arises 
naturally in various physical contexts, such as the diffusion of particles, slow relaxations in glasses, 
and scalar phonon localization. Using a combination of a renormalization group procedure and a 
direct moment calculation, we find the eigenvalue distribution density (i.e., the spectrum), for low 
densities, and the localization properties of the eigenmodes, for arbitrary dimension. Finally, we 
discuss the physical implications of the results. 
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Application of the theory of random matrices whose 
elements are independent Gaussian variables has proven 
to be rich mathematically and relevant for many physi- 
cal systems In this Letter we study a different class 
of random matrices where the «, j'th element is a func- 
tion of the Euclidian distance between pairs of points 
whose positions are chosen randomly and uniformly in a 
d-dimensional space. It is natural that in cases where the 
matrix element is related to an overlap between localized 
quantum-mechanical wavefunctions, the dependence on 
the distance will be exponential, i.e., Ay = e~ rij ^, with 
£ being the localization length [2j. 

The exponential matrix is an appropriate model for 
various physical systems, in this Letter we will concen- 
trate on its application to glasses relaxing to equilibrium, 
a particle diffusing in random environment and localiza- 
tion of phonons. Most of the results are derived at the 
low density limit, when e = £/r„„ <C 1, with r nn being 
the average nearest neighbor distance. To understand the 
properties of these systems one need to find out the dis- 
tribution density P(A) of the eigenvalues A as well as the 
structure of the eigenmodes. An intuitive picture of the 
problem arises in the application to phonon localization 
with springs constants i£y that depend exponentially on 
the Euclidean distances between the masses; we therefore 
use the phonon terminology: eigenmode. 

The low density limit allows us to find P(A) analyti- 
cally employing a direct calculation of its moments, see 
Eq. (J2J and the Supplementary Material (SM). We find 
that -P(A) ~ 1/A in all dimensions over a broad range of 
A's. While in one dimension the normalization of P(A) is 
assured by an intcgrable power-law divergence at eigen- 
values close to zero, for higher dimensions there is a peak 
related to a finite cutoff, cf Fig. [T] We use a logarith- 
mic scale to plot P[log(— A/2)] in order emphasize the 
deviations from the 1/A distribution. 

To comprehend the structure of the eigenmodes we use 
a renormalization group (RG) approach for random sys- 
tems that was developed in the context of spin chains 



At each RG step, we choose the stiffest spring. 
Since the spring is large by construction, after finding 
the eigenvalue associated with the stiffest spring we can 
glue together the two masses at its ends creating a larger 
mass. At the next RG step we choose again the stiffest 
spring among those who remain. In this way the large 
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FIG. 1: Comparison between the exact numerical diagonal- 
ization (circles), and different theoretical approaches: The 
stars (green) show the results of a renormalization group ap- 
proach, where, if one considers the problem as finding the 
normal modes of a spring network, the springs and masses 
are renormalized. The solid line (red) depicts the analytical 
results of Eq. ©, and Eq. (11) of the SM. The numerical 
results shown are for 1, 2 and 3 dimensions, with N=1000, 
averaged over 1000 realizations. The points were chosen in a 
line, square or box of side one, and e = 0.1, corresponding to 
£ = 10" 4 in ID, £ = 0.0032 in 2D and £ = 0.01 in 3D. Notice 
that the graph shows the distribution density of the logarithm 
of the eigenvalue, which eliminates the governing 1/A depen- 
dence, and allows us to observe clearly the deviations from it. 
To demonstrate this, the horizontal dashed line (blue) shows 
an exact 1/A distribution. No fitting parameters are used. 
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cigcnmodcs are built initially by pair of masses, but as 
the RG process progress larger clusters of masses form 
eigenmodes with smaller eigenvalues. This behavior is 
demonstrated in Fig. [5] 

Since in the RG step a spring is chosen regardless of the 
masses, the flow of the mass distribution is independent 
of that of the springs. Namely, we will get statistically 
the same distribution of clusters as a function of the RG 
step by a coagulation process where we randomly choose 
two masses and glue them together at each step. This 
observation allows us to find the mass distribution as a 
function of the number of RG steps (see Eqs. (7) and 
(8) of the SM). Together with the exact result for the 
eigenvalue distribution, we can find the dependence of 
the cluster mass on the eigenvalue n c (A), cf Eq. ([5]). 

We will now present in more details the properties of 
the exponential random matrix model, its applications, 
and the derivation of P(X) and n c (A) 

Model and relevant physical problems.- N points are 
chosen randomly in a d dimensional cube, and is de- 
fined as the Euclidean distance between points i and j. 
We define a matrix A, as follows: Aij = e^ rij ^, for i ^ j 
and An = — ^2j^i Aij , the latter definition expressing 
a conservation law in the physical problem [7H9(. We 
shall be interested in determining P(X), the probability 
density of eigenvalues of the matrix A, for low densities 
(small values of e). Since the matrix is hermitian, it is 
clear that all its eigenvalues are real, and it can also be 
proven that they are negative [1, [l(| . 
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FIG. 2: Demonstration of the structure of the eigenmodes. 
N = 5000 points were chosen randomly and uniformly in a 
two dimensional box, with e = £/r nn = 0.1. The eigenmodes 
become more delocalized as the eigenvalues approach zero, a 
condition made quantitative in Eq. ©, through the use of 
the RG approach. Eigenmodes are comprised of clusters of 
points, localized in real space. 



This model is relevant for problems from various fields 
of physics, but for now, we choose to focus on scalar 
phonon localization. It will be useful to have this problem 
in mind when we discuss the RG calculation. 

When studying normal modes of a collection of equal 
masses m connected by harmonic springs, one has to find 
the eigenmodes of a matrix A, where {^4^ } are the spring 
constants, and due to momentum conservation the sum 
of columns vanishes. The eigenvalues are related to the 
frequencies by mu> 2 = —A, where we can choose m to be 
unity, for convenience. The above model can be used to 
study phonons in a disordered lattice, where due to the 
randomness in the matrix elements (i.e., a distribution of 
spring constants), the oscillating modes may turn out to 
be localized in sj 
phonon model [1! 
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Notice that this is a 'scalar' 
| , where the vectorial properties of the 
modes have not been taken into consideration (and thus 
the matrix has N and not Nd modes). 

For phonons, the density-of-states of the vibrations 
P(u>) = 2mujP(X = —muj 2 ) determines the thermal prop- 
erties of the system [l3| . It is interesting to note that our 
result for P(A) implies P(oS) ~ l/u> up to corrections de- 
pending on the dimensionality 

Derivation of the eigenvalue distribution.- The expo- 
nential matrices model in one dimension has peculiar 
spectral properties [TUtHI- The gi yes a simple, non- 
rigorous derivation for the low density limit of the eigen- 
value distribution in one-dimension, which turns out to 
be a power-law. The argument also shows that in one 
dimension the modes are localized, with typical size de- 
pending on A as a power-law. 

We will now develop a general formula for the eigen- 
value spectrum at arbitrary dimension, in the low density 
limit. We shall see that in one dimension we obtain a 
power-law, while in higher dimension we get logarithmic 
corrections to a distribution P(X) ~ 1/A . 



Using the following sum: 



J P(X)X k dX 



E 



Ai 1 Ai 2 i 3 ...Ai k i x 



(1) 

one can calculate the fc'th moment of the probability 
density. In the SM we find that in the low density 
approximation e — ► 0, Ik ~ — (— 2) k ~ 1 d\Vd(e/k) d , with 
Vd = Ti d / 2 /T(d/2 + 1) the volume of a d dimensional unit 
sphere. 

Since a distribution is fully determined by its moments 
under certain conditions which are fulfilled in our case 
[l6| . it suffices to find a distribution which yields these 
moments. It can be directly checked by performing the 
integrals that the following probability density does this: 



P(A) 



e d dy d /2[-log(-A/2)] d - 1 c-^ ied (- lo s(- A /2)) d 



(2) 
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The cumulative C(A) of this distribution takes a par- 
ticularly simple form: 



C(A) = J P(X)d\ = e 



- p --^ d (-log(-A/2)) d 



(3) 



In one-dimension we find a power-law divergence in 
the distribution density, while in a dimension d > 1 we 
see from Eq. ([2]) that there is a finite cutoff at small 
eigenvalues. Fig. [T]comparcs this formula with numerical 
diagonalization, for the one, two and three dimensional 
cases, as well as the RG procedure. 

It is important to distinguish the above result from, 
for example, those of Ref. [H, which relates the eigen- 
modes to isolated pairs of points. There, an uncontrolled 
approximation was used, connecting pairs of nearest- 
neighbors (a procedure which is not always well-defined) . 
While giving the correct qualitative dependence, it dif- 
fers from Eq. ^ by a factor of 2 in the exponential, 
as well as in the normalization factor. As will be shown 
via the renormalization group approach in the next sec- 
tion, this difference reflects important differences in the 
underlying physics, and in the structure of the eigen- 
modes: the high frequency cigenmodes indeed consist of 
pairs of nearest-neighbors masses, and for them the two 
approaches coincide. However, the low-lying modes are 
comprised of an increasing number of masses, diverging 
as the frequency approaches zero. For any finite value of 
e, Eq. ([2]) might be violated for sufficiently small A, while 
keeping the corrections to all moments of the distribution 
sufficiently small. 

Renormalization group approach.- Following let 
us consider the renormalization group (RG) approach to 
this problem. This will enable us to understand the lo- 
calization properties of the eigenvectors. We recall the 
scalar phonon localization picture, where the points rep- 
resent masses and the matrix-elements represent spring 
constants. At each RG step, we choose the largest spring 
K: due to the broadness of the spring distribution, it will 
not be affected much by the neighboring springs, and 
therefore it will contribute a frequency oj 2 = — = —A, 
where a = mi ''" 2 is the reduced mass of the two masses 

f rai+mj 

it connects, mi and ui2- Initially, all masses are equal. 
Notice that the mechanical intuition tells us we should 
choose the largest spring at each step, and not necessar- 
ily the largest frequency, i.e., the choice of the springs 
does not depend on the masses. Since this spring is large 
by construction, we can now 'glue' the two masses to- 
gether, and create a single mass mi + m.2- We also have 
to renormalize the springs attaching the new mass and 
all other masses: a reasonable way to do so is to take 
each of the springs between a mass m and the new (joint 
mass) as the sum of the two springs between the mass 
m and each of the masses mi and mi- Clearly we will 
obtain smaller frequencies and springs and larger masses 
as the RG progresses. 



As is shown in Fig. (fTJ), this simple scheme cap- 
tures the essential physics, with no fitting parameters. 
As mentioned, the reason the method works so well is 
the broadness of the 'spring' distribution: for a one 
dimensional case, for example, the distribution of the 
nearest-neighbor spring constants (which can be calcu- 
lated directly from the exponential distribution of the 
distance intervals) follows a P(K) ~ l/if 1-6 , where 
e — > in the low density limit. Notice that for the 
one-dimensional case the RG procedure would choose 
exactly these nearest-neighbor springs, by construction. 
Thus, were we to neglect the mass RG we would obtain 
P(X) ~ 1/A 1_e , which recovers the low density result 
mentioned earlier. In the SM we show how one can in- 
corporate the mass RG to correct the one-dimensional 
result also for higher densities. 

Since at each RG step a spring is chosen regardless 
of the masses, the flow of the mass distribution is com- 
pletely independent of that of the springs. As mentioned, 
the results of the mass RG are applicable for any di- 
mension. We will now show that using the probability 
density we obtained, Eq. (fj)), with the results of the 
mass RG, we can understand the localization properties 
of the eigenvectors. The SM shows how we can find the 
distribution of masses at a given stage of the RG pro- 
cess, which turns out to be approximately an exponen- 
tial distribution e~ m ^ m \ with (m) changing as the RG 
process evolves, corresponding to the formation of larger 
and larger clusters. This implies there is a typical mass 
(m) at each instance. Since at each step of the RG pro- 
cess the number of clusters decreases by one, after k steps 
the average cluster mass is given by: (to) = N/(N — k), 
which is also the typical size of the cluster, n c . On the 
other hand we can find the relation between the RG step 
k and the eigenvalue A: In the RG process at each step 
the number of masses is decreased by one, and the corre- 
sponding cigenmode recorded, starting with the highest 
springs (and eigenvalues). At the stage of the RG flow 
corresponding to an eigenvalue A, the number of masses 
left, N — k, can be calculated using Eq. ([3]): 



^P{X)dX = N(l - C(A)). 



(4) 



Combining this equation with that for n c = (to), we 
find that n c depends on the eigenvalue as: 



n c (X) ~ 1/C(X) 



e 2 



i (-log(-A/2)) d 



(5) 



As we go to zero eigenvalue, the size of the eigenmodes 
diverges, as demonstrated in Fig. ^ on a particular 
example. Notice that for the case d = I we recover 
the power-law relation between localization length and 
A mentioned earlier, and related to Refs. 



14) and [I5|. 



Eq. (j4]) also allows us to count Af(n c ), the number of all 
clusters containing more than n c masses: their number is 
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N J?P(X)dX = NC(\), where we know the dependence 
of A on n c through Eq. (0. This gives: Af(n c )/N = l/n c . 

Physical implications.- The mathematical model pre- 
sented is relevant also for other physical problems besides 
phonon localization discussed above. For example, one 
can consider the hopping of a particle in a random envi- 
ronment, where Aij describes the transition probability 
of a particle from site i to site j. If we define the prob- 
abilities of the particle to be at the different sites by a 
vector p, then -§ = Ap. In disordered systems, the tran- 
sition rate often depends exponentially on the distance 
[I?}- Another physical example of relevance, is the study 
of relaxations in glasses. Under certain approximations, 
this can be mapped to the study of cigenmodes of a class 
of random matrices related to the one described above 
[l8j and in fact this has been the original motivation for 
this study. In both the cases of the diffusion problem 
and the relaxations in glasses the Laplace transform of 
the distribution density plays an important role: in the 
hoppingproblem, it gives the probability to remain in the 
origin [19(, while for the glass relaxation it corresponds 
to the time dependence of the relaxation [H, [2(| ■ 

Upon taking the Laplace transform P(t) of the distri- 
bution density in one-dimension (with argument i), we 
obtain P(t) ~ t~~ . This implies that the diffusion in 
this case is anomalous [l7, 21, 22|: the probability to re- 
main in the origin does not decrease as 1/ s/t, as is the 
case for normal diffusion, but as a smaller power (the 
particle tends to be more localized). If we assume that 
after a time t the particle has spread over a distance r, 
it is reasonable to assume that p(t)r ~ 1, showing that 
r ~ t't+t . In the case e — > 0, this can be approximated 
as P(t) ~ C — log(t). This type of behavior has been 
experimentally observed in various types of glassy sys- 
tems (23M25| . In higher dimensions, one obtains in the 
low density limit P(t) = 0[log(t)], with <j) a polynomial 
of degree d. This will be elaborated on in future works. 

Summary.- We presented here a model of random ma- 
trices which captures the interesting physics of various 
different systems. After introducing the model, we found 
the eigenvalue distribution density in the low-density 
limit (Eq. ((2J) and the localization properties of the 
eigenmodes (Eq. (|5|l) using a direct moment calculation 
as well as a renormalization group approach. Our results 
for the spectrum agree with the ID case LJ, |l5j, and 
with the exact numerical diagonalization. While in one 
dimension it is known that for an infinite system there 
will be an (integrable) power-law divergence of the spec- 
trum, in a higher dimension d we found that there is a 
finite cutoff, where e = £/r nn is the small parameter of 
the theory. We used the RG approach to show that the 
eigenmodes are localized, and to find a relation between 
the spatial extent of an eigenmode and the corresponding 
eigenvalue, implying that this size diverges as the eigen- 
value approaches zero. 

We discussed the application of the model for various 



physical problems, such as relaxations in glasses, diffu- 
sion of particles in random media, and localization. In 
the future, it would be fascinating to understand also the 
crossover or phase-transition to the high density or low 
disorder regime, which should present different physics. 

We thank M. Aizenman, D. Cohen, M. E. Fisher, O. 
Hirschbcrg, B. Nedler, W. Schumacher, V. Vitclli and 
O. Zeitouni for important discussions. This work was 
supported by a BMBF DIP grant as well as by ISF and 
BSF grants and the Center of Excellence Program. 



SUPPLEMENTARY MATERIAL 



1. A NONRIGOROUS DERIVATION IN 
ONE-DIMENSION 



It is known that for an isolated system the center-of- 
mass motion gives a zero eigenvalue mode, and thus in 
our model a cluster of points largely separated from the 
rest of the points will give an exponentially vanishing 
eigenvalue. Let us consider the statistics for finding such 
a cluster, isolated on both sides by a distance r (here 
we explicitly use the one dimensional character of the 
problem) . 

By perturbation theory the vanishing eigenvalue will 
be shifted by a value of order |A| ~ e~ r ^, where £ was 
the decay length of the matrix elements. We choose the 
points from a uniform distribution, hence the distances 
between them are distributed exponentially. Therefore, 
if we move (e.g. to the right) along the one dimen- 
sional chain, the probability of encountering a gap of 
size f larger than r is P(r > r) = e~ r / r ™ n with r nn 
being the average distance between the points. If we 
proceed to move along the chain we are certain to en- 
counter a gap larger than r eventually. Thus, the proba- 
bility density of finding a cluster isolated from both sides 
by a distance r (obtained by differentiating the cumula- 
tive distribution) is p(r) = e~ r / r ""/r ml . Using the re- 
lations, P(X)d\ = p(r)dr we find P(X) oc l/A 1 ^, with 
e = £/r„„ being a small parameter if the density of points 
is low enough. This analysis also allows us to find the 
length of the cluster, since the typical number of points 
we have to pass until we encounter the second gap is 
l(X) = 1/p oc e r l Tnn ~ A -1 / c . This reproduces the low 
density limit of a similar problem [l4|, EH , albeit with 
nearest-neighbor matrix elements only (indeed, in the low 
density limit, in one dimension, we do not expect next- 
nearest neighbor couplings to modify the results). The 
relation l(X) that was obtained implies that the spatial 
extent of the modes diverges as the eigenvalue approaches 
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2. CALCULATION OF THE MOMENTS OF THE 
EIGENVALUE PROBABILITY DENSITY 

We would like to calculate the fc'th moment of the 
probability density using the following sum: 

Jp(X)X k dX = ±( ]T A^ 2 A l2 ^...A lk , n ) (6) 

ii,i 2 ..ik 

To gain intuition, let us consider the first moment: 
/ P(X)XdX = Ai.i) ■ By construction, the magni- 

tude of the diagonal elements is the sum of all the other 
elements in the columns, i.e., sum of the hopping matrix 
elements to all neighbors, close and far. To the lowest or- 
der in e, we can consider only the nearest-neighbor of each 
point, as we shall shortly demonstrate. Under this ap- 
proximation the first moment can be readily calculated: 
h = J P(\)\d\ ps - J™ P nn {r)e- r /Sdr, where P nn (r) is 
the nearest-neighbor distance distribution, which is given 
by P nn (r) = f^(r/r nn ) d - x e- v ^l r ^\ where V d is the 
volume of a d dimensional sphere and r nn is the average 
nearest-neighbor distance. In the low density limit this 
gives I\ ps — d\Vd€ d - Calculating the corrections to this 
formula due to all neighbors, one obtains a correction of 
the order of 0(e 2d ), which is indeed of higher order in the 
small parameter e. The crux of the matter here is that 
due to the exponential function, the contributions to the 
integral arise mainly from rare pairs which are very close 
to each other, with a distance of order £ between them. 
For such pairs the probability that a third point will be 
in the vicinity of the two points is negligible. 

For the fc'th moment the general term contributing will 
be a combination of hops to the nearest-neighbor and 
loops that stay in either the original site or the nearest- 
neighbor, as illustrated schematically in Fig. G3 Thus, 
there would be 2 fe ~ 1 contributing diagrams, each of equal 
contribution. Therefore the fc'th moment will be given by 
h = JP(X)X k dX ps -(-2) fc - 1 / O °^n„(r-)e- fcl -/^ r . This 

yields in the low density limit: I k ps -(-2) k - 1 d\V d (e/k) d . 
These moments correspond to the probability density 
given in Eq. (2), as can be verified by a direct calcu- 
lation of its moments. 




FIG. 3: An example of a diagram contributing to the fc'th 
moment in the low density limit. 

It is interesting to note that although the leading or- 
der contribution to the fc'th moment was made by a 
small subset of configurations, close pairs, the resulting 



distribution captures the probability density of eigenval- 
ues whose corresponding eigenmodes are actually clus- 
ters. The clusters can understood by the renormalization 
group approach, discussed in the manuscript. It is gener- 
ally true that for large fc's, the moments are dominated 
by the largest (i.e., most negative) eigenvalues, which in 
our case are exactly the close pairs (which are nearly 
eigenmodes). The above calculation shows that in the 
low density limit also the low moments (e.g: I\, above), 
are dominated by the close pairs. Thus, we can find all 
moments and the complete distribution, by studying the 
properties of the close pairs, in a similar fashion to an 
analytical continuation, by which studying the function 
in a small region determines its behavior further away. 

In the case of matrix elements decaying not exponen- 
tially but as some function f(r) of the distance, we can 
repeat the previous argumentation. If f(r) decays, for 
example, faster than exponential, than surely the contri- 
bution of the diagrams neglected will be smaller still, thus 
justifying the method. In this case the fc'th moment will 
be: I k = jP(X)X k dX ps -(-2)*" 1 / °° P nn (r)f(r) k dr. It 
can be checked that the cumulative of the associated dis- 
tribution in this case is given by: 

C(X) = J°P(X)dX = e '^t arl{ ' X/2))d , (7) 

where / _1 is the inverse function of f(r). As an exam- 
ple, applying this formula for gaussian matrix elements 
in two-dimensions gives a power-law distribution density, 
P(X) ~ I/A 1 -^ 2 / 2 . 

3. SOLUTION FOR THE RENORMALIZATION 
OF THE MASSES 

As explained in the manuscript, the RG procedure 
chooses the largest springs at every instance. With re- 
gard to the masses this is a random process not depend- 
ing explicitly on the masses, as long as we neglect effects 
arising from the larger surface area of the clusters. As 
is shown in the manuscript, these do not affect (to). We 
shall now analyze the resulting flow of the mass distribu- 
tion, under this approximation. Initially, we have a set 
of unit masses, and at each step we choose two masses 
at random and combine them. If we denote the average 
number of masses of mass m in the fc'th stage as Nf^ (at 
which step we have N — k masses left), we therefore have: 

N m -N m - w -^+ ^ _____ ( 8 ) 

m'=l..(rn— 1) 

Indeed, the probability that a chosen mass is of mass 

m i is M-k ■ This relates to the first term, which accounts 
for the average reduction in the number of masses of mass 



6 



m (the factor of two arising from the fact that we choose 
two masses at each stage) . The second term accounts for 
all the possible ways that two smaller masses can combine 
to contribute to the number of masses of mass m, where 
again the same form for the probabilities enters. 

In order to move from numbers to probabilities, we 
have to divide by the current number of masses, which 
is N — k — 1 (since at each RG step, when two masses 
are combined to create one mass, we lose one mass). It 
is useful to write, for N — k 3> 1: 



N-k-1 N - k (N - k) 2 ' w 

Denoting by the probability to have a mass of size 

m after k springs are eliminated, we have P^ n = 3", , 
leading to: 

pk o pk 1 

^ pk | rn m i \ pk pk 

m ~ m^ N _ k N _ k ^ N _ k m' m- 

m' — l..(m — 1) 

(10) 

Thus, we finally obtain: 

AP m W -Jj— £ + £ Pm'Pm-m'i (H) 

m'=l..(m— 1) 

where the second term of Eq. (|9]) changed the factor 
of two appearing in Eq. ((5J) to be one. 

It is now straightforward to take the continuous limit, 
where it is convenient to define a fictitious time related 
to the step number as t — — log . This leads to the 
following intcgro-diffcrcntial equation for the time evolu- 
tion of the masses: 



the result of the moment calculation. Eq. (TT3]) shows 
that the mass distribution at a given time is given 
by V(m,t) = e~ e m e~ t . From this, we can calcu- 
late (I/to) = j v{ - m ^ dm. Approximating the expo- 
nential as an effective cutoff at m c = e*, we find that 
(1/m) ~ log(m c )c _t ~ te _t , where we took a lower cut- 
off of order unity for the mass. 

Using the relation between t and k we find that: 

(1/m) ~ (N- k)/N, (14) 

where N — k is the number of masses left after k RG 
steps. 

For the one dimensional case we can combine this re- 
sult with the previously obtained renormalization of the 
springs, to obtain the mass correction. Since the eigen- 
modes obey w 2 = Kj with \i the reduced mass of the 
two current masses, the springs eliminated at this point 
of the RG flow, whose distribution was shown to obey 
P{K) ~ l/K 1 - € , contribute on average to 

- A = u? = K{l/m! + l/m 2 ) = 2K{l/m) = 2K 1+e . 

(15) 

Calculating the distribution of A arising from the known 
spring distribution P{K) gives 

P(A)~1/Atri, (16) 

where the lower cutoff is still A = but the higher cutoff 
is A = —2. Thus, the effect of the mass RG in one di- 
mension is changing P(A) ~ 1/A 1_e to I/A1+ 7 , which is 
negligible for small e. For low densities the RG approach 
and the derivation via the calculation of the moments 
give identical results, as they should. 



8V(m t) r n 

— = -P(m,t) + / V(m',t)V(m-m',t)dm', 
ot Jo 

(12) 

where the time t is related to step number k by the equa- 
tion t = — log[(N — k)/N]. This is a particular case of 
the Smoluchowski coagulation equation p6| . It can be 
verified that: 



V{m,t) = e 



(13) 



is a solution. We expect this solution to be relevant 
for t ^> 1, when the slightl y d ifferent initial conditions 
it obeys are not important [271 ] (initially, all masses are 
equal, not exponentially distributed). 



4. RG SOLUTION IN ONE DIMENSION 

In the following we show how in one dimension one 
can find P(A) with the RG approach, without using 
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